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The  kinetic  parameters  at  end-of-life  of  a  material  test  reactor  fuelled  with  low  enriched  uranium  fuel 
were  calculated.  The  reactor  used  for  the  study  was  the  IAEA’s  10  MW  benchmark  reactor.  Simulations 
were  carried  out  to  calculate  core  excess  reactivity,  neutron  flux  spectrum,  prompt  neutron  generation 
time  and  effective  delayed  neutron  fraction.  Nuclear  reactor  analysis  codes  including  WIMS-D4  and  CITA¬ 
TION  were  employed  to  carry  out  these  calculations.  It  was  observed  that  in  comparison  with  the  begin- 
ning-of-life  values,  at  end-of-life,  the  neutron  flux  increased  throughout  the  core,  the  prompt  neutron 
generation  time  increased  by  3.68%  while  the  effective  delayed  neutron  fraction  decreased  by  0.35%. 
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1.  Introduction 

A  large  number  of  research  reactors  around  the  world,  which 
were  fuelled  with  HEU  based  fuels  having  uranium  enrichment 
around  90%  in  235U  isotope  have  been  converted  to  use  LEU  based 
fuels  having  uranium  enrichment  of  20%  in  235U  isotope,  since  20% 
fuel  enrichment  is  an  isotopic  barrier  for  weapon  usability  (Glaser, 
2005).  The  IAEA  has  devised  a  standard  benchmark  MTR  (IAEA- 
TECDOC-233, 1980)  in  order  to  facilitate  reactor  conversion.  Many 
theoretical  calculations  have  been  performed  and  reported  in 
various  documents  (IAEA-TECDOC-233,  1980;  1AEA-TECDOC-643, 
1992). 

All  the  calculations  reported  in  the  IAEA’s  guidebooks  (IAEA- 
TECDOC-233,  1980;  IAEA-TECDOC-643, 1992)  deal  with  the  reactor 
behaviour  at  the  beginning  of  reactor  core  life.  However,  the  reac¬ 
tor  remains  at  the  beginning-of-life  (BOL)  just  for  a  while  and  as 
soon  as  the  fission  reaction  starts,  the  isotopic  concentration  in 
the  fuel  changes.  Since  most  of  the  neutronic  parameters  depend 
on  the  fuel  composition,  these  parameters  also  change  with  the 
change  in  the  fuel  isotopic  concentration.  Hence,  need  is  always 
felt  to  find  different  reactor  parameters  at  the  end-of-life  (EOL). 

The  work  presented  in  this  paper  deals  with  the  calculation  of 
the  neutron  flux,  prompt  neutron  generation  time  and  effective 
delayed  neutron  fraction  the  10  MW  IAEA  benchmark  reactor 
(IAEA-TECDOC-643,  1992)  using  LEU  fuel  at  the  EOL. 
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2.  Reactor  description 

The  reactor  analyzed  is  the  same  one  utilized  for  the  benchmark 
problem  solved  in  IAEA-TECDOC-233,  with  the  water  in  the  central 
flux  trap  replaced  with  a  7.7  cm  x  8.1  cm  block  of  aluminum  con¬ 
taining  a  square  hole  5.0  cm  on  each  side  (IAEA-TECDOC-643, 
1992).  Description  of  the  low  enriched  uranium  core  of  the  reactor 
is  given  in  Table  1.  The  core  configuration  and  burn  up  of  fuel  ele¬ 
ments  in  percentage  of  loss  of  the  number  of  initial  235U  atoms  at 
BOL  is  given  in  Fig.  1  while  that  of  EOL  is  given  in  Fig.  2.  Other  de¬ 
tails  can  be  found  in  the  reference  documents  (IAEA-TECDOC-233, 
1980;  IAEA-TECDOC-643,  1992). 

3.  Analysis  procedure 

3.1.  Reactor  simulation  codes 

The  WIMS-D4  (Halsall,  1980)  code  was  used  for  the  genera¬ 
tion  of  group  constants  for  various  core  regions  while  CITATION 
(Fowler  et  al„  1971)  was  used  to  perform  global  core  calculations. 
Detailed  description  of  these  codes  can  be  found  in  related 
material. 

3.2.  Simulation  methodology 

The  CITATION  code  was  used  to  calculate  various  core  parame¬ 
ters  like  ke ff,  neutron  fluxes  and  adjoint  fluxes.  The  core  was 
simulated  in  the  x-y-z  geometry.  All  control  rods  were  assumed 
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Table  1 

Data  for  low  enriched  uranium  core  of  IAEA  10  MW  benchmark  reactor  (IAEA- 
TECDOC-643,  1992). 

Parameter  Value 


Active  core  height 

Extrapolation  length 

Space  at  the  grid  plate  per  fuel 
element 

Fuel  element  cross-section 

Thickness  of  support  plate 

Number  of  fuel  plates  per  fuel 
element 

Number  of  fuel  plates  per  control 
element 

Identification  of  the  remaining 
plate  positions  of  the  control 
element 


of  the  LEU  (UAljc-Al) 


Xenon-State 

Nominal  coolant  flow  rate  (m3/h) 
Coolant  inlet  temperature  (°C) 
Hot  channel  factors 


60.0  cm 

the  cosine-shaped  flux  goes  to  zero) 

7.6  cm  x  8.05  cm  including  support 
plate 

6.3  cm  X  0.051  cm  x  60.0  cm 
0.475  cm 

23  identical  plates,  each  0.127  cm  thick 
17  identical  plates,  each  0.127  cm  thick 

Four  plates  of  aluminum,  each  0.127  cm 
thick  in  the  position  of  the  first,  the 
third,  the  twenty-first,  and  the  twenty- 
third  standard  plate  position;  water  gaps 
between  the  two  sets  of  aluminum 
plates 

a. Enrichment  90  w/o  U-235  in  HEU,  20 
w/o  U-235  in  LEU 

b. 390  g  U-235  (ptJ  =  4.40  g/cm3)  per  fuel 
element  (23  plates) 

c.  72  w/o  of  uranium  in  the  UA1X-A1 

d. Only  U-235  and  U-238  in  the  fresh  fuel 
10MWth 

Homogeneous  Xenon  content 
corresponding  to  average-power-density 

38 

a. Radial  x  local  power  peaking 
factor = 1.4 

b. Axial  power  peaking  factor*  1.5 

c. Engineering  factor  =  1.2 

d. Overpower  factor  =  1 .2 


to  be  fully  withdrawn.  The  fuelled  and  non-fuelled  portions  of  each 
standard  and  control  fuel  element  were  modeled  separately.  The 
WIMS-D4  code  was  used  for  computation  of  macroscopic  absorp¬ 
tion  cross-section  (Z„),  the  v-fission  cross-section  (vZ/),  the  diffu¬ 
sion  coefficient  (D),  the  scattering  matrix  (ZSig  ,g>)  and  the  fission 
spectrum  for  all  groups.  These  data  are  required  by  CITATION  as  in¬ 
put.  Five  energy  groups  were  used  to  perform  multi-group  calcula¬ 
tions.  These  groups  are  given  in  Table  2.  The  simulation 
methodology  has  been  fully  described  and  validated  in  our  recent 
works  (Muhammad  and  Majid,  2008,  2009). 


Outside  boundary 


Table  2 

Energy  groups  used  for  macroscopic  cross-section  generation  by  WIMS-D4. 

Group  no.  E0  (eV)  EL  (eV)  Group  type3 *  Flux  typeb 

1  10.0  X  106  0.821  X  106  Fast  Fast 

2  0.821  X  106 *  5.530  x  103  Resolved  resonances 

3  5.530  x  103  0.625  Unresolved  resonances  Epithermal 

4  0.625  0.14  Thermal  Thermal 

5  0.14  0  Thermal 


3  Muhammad  and  Majid,  2008. 
b  IAEA-TECDOC-233,  1980. 


4.  Results  and  discussions 

4.1.  The  keff  at  EOL 

The  CITATION  code  gives  the  value  of  keJj- directly.  This  value  for 
the  given  core  configuration  at  EOL  was  found  to  be  0.999560 
showing  that  the  reactor  fuel  had  actually  come  to  the  end  of  its 
life  (Table  3). 


4.2.  Neutron  flux  spectrum 

The  shape  of  neutron  fluxes  of  all  the  groups  is  same  (Fig.  3)  at 
BOL  and  EOL  but  the  magnitude  of  flux  of  all  the  three  groups  in¬ 
creases  at  EOL  (Fig.  4).  The  thermal  neutron  flux  at  the  mid  of 
the  central  flux  trap  also  increases  (Table  4)  to  2.2316  x  1014/ 
cm2  s  from  2.1780  x  1014/cm2  s  showing  an  increase  of  2.48%  from 
the  BOL  value  calculated  using  same  method  and  technique 
(Muhammad  and  Majid,  2008).  The  reason  behind  this  increase  is 
that  the  quantity  of  fuel  decreases  with  burn  up  decreasing  the  fuel 


Table  3 

Comparison  of  different  parameters  at  BOL  and  EOL. 

Parameter  BOL  EOLf 

Ref.  I3  Ref.  2b  Ref.  3C  Ref.  4d  Ref.  5e 


keff  1.01823  -  -  1.02462  1.018273  0.999560 

A  (ps)  43.74  -  44.39  45.42  44.03  45.65 

jleg  0.007275  0.00732  0.007219  0.007351  0.007185  0.007160 


3  IAEA-TECDOC-643,  1992  (Appendix  G-l). 

b  IAEA-TECDOC-643,  1992  (Appendix  G-2). 

c  IAEA-TECDOC-643,  1992  (Appendix  G-3). 

d  IAEA-TECDOC-643,  1992  (Appendix  G-4). 

'  Muhammad  and  Majid,  2008. 

f  This  work. 
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The  value  of  A  at  EOL  was  found  to  be  45.65  ps  (Table  3)  show¬ 
ing  an  increase  of  3.68%  over  its  value  at  BOL  (Muhammad  and 
Majid.  2008). 

During  the  reactor  operation,  the  total  amount  of  fissile  isotopes 
(235U  and  239Pu)  decreases  linearly  with  burnup  (Fig.  5)  (The  con¬ 
centration  of  241  Pu  is  too  low  to  have  any  significant  effect).  This 
linear  decrease  is  seen  in  the  value  of  v  E/  as  shown  in  Fig.  6 
(especially  in  groups  4  and  5  where  most  of  the  fission  takes  place), 
which  also  decreases  linearly  with  fuel  burnup.  The  increased  va¬ 
lue  of  A  reflects  the  decreased  value  of  v  E  /  and  increased  ther- 
malization  of  the  neutron  flux  at  EOL. 

4.4.  The  effective  delayed  neutron  fraction  fteff 


0  10  20  30  40  50 

x-Distance  (cm) 

Fig.  4.  Ratio  of  neutron  fluxes  at  EOL  to  that  at  BOL. 


Table  4 

Comparison  of  neutron  flux  at  mid  of  central  flux  trap  at  BOL  and  EOL  (1014  n/cm2  s). 
Flux  type  BOLa  EOLb  %  increase  from  BOL  to  EOL 

Fast  1.789  1.806  0.95 

Epithermal  0.963  0.975  1.25 

Thermal  2.178  2.232  2.48 


a  Muhammad  and  Majid,  2008. 


number  density  and  since  the  power  remains  constant  at  10  MW, 
the  neutron  flux  increases. 

4.3.  The  prompt  neutron  generation  time  A 

The  prompt  neutron  generation  time  {A)  was  calculated  from 
the  following  relation  (Ott  and  Neuhold,  1985). 


where  (1)  is  the  average  inverse  neutron  velocity  given  by  (Ott  and 
Neuhold,  1985): 

m  L  ,4  ^Pir,E)dEdv  siU&H  ,,, 

W  fv  JMr  E)dEdV  Eg=t<?>g  U 

where  <pg  is  the  group  flux  and  Vg  is  the  group  average  neutron 
velocity,  and 

v  /vJe  v5>(r,E)  (p(r,E)dEdV  Eg.,vE/g<Pg 
jv  fE  <p(r.  E)dEdV  "  Eg=i<Pg 


The  peff  was  calculated  from  the  following  relation  (IAEA-TEC- 
DOC-643,  1992). 


s  eUva  ’ 


(4) 


where  i  is  an  index  for  fissionable  isotopes,  yg  is  g-group  fraction  in 
typical  delayed  neutron  spectrum,  Xg  is  g-group  fraction  of  fission 
spectrum,  /('  is  i-nuclide  delayed  neutron  fraction,  <p+  is  core-aver¬ 
age  adjoint  flux  for  group  g,cp'g  is  core-integrated  g-group  ordinary 
flux,  (vEw>)t  is  total  source  of  fission  neutrons.  The  values  of  y’g 
and  [i1  were  taken  from  reference  (Keepin,  1965.)  while  the  remain¬ 
ing  required  data  were  obtained  from  WIMS  and  CITATION.  The  va¬ 
lue  of  at  EOL  was  found  to  be  0.007160  (Table  3)  showing  a 
decrease  of  0.35%  over  its  value  at  BOL  (Muhammad  and  Majid, 
2008).  The  amount  of  different  isotopes  changes  from  BOL  to  EOL 
(Fig.  5)  resulting  in  change  in  the  value  of  peff.  Since  the  amount 
of  Pu  isotopes  and  their  change  in  concentration  with  burnup  is 
small  as  compared  to  the  amount  of  235U  isotope,  the  change  in 
the  value  of  is  also  small. 


2.5E-03 
2.0E-03 
1.5E-03 
1.0E-03 
5.0E-04 

0.0E+00  ..... 

0  10  20  30  40  50 


Burnup  (%  of  initial  235U  atoms) 


Fig.  5.  No.  densities  of  different  isotopes  with  fuel  burnup. 


the  required  data 


obtained  from  WIMS  and  CITATION. 


Fig.  6.  Variation  of  v  2/  with  fuel  burnup. 
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5.  Conclusions 

The  kinetic  parameters  of  an  MTR  fuelled  with  LEU  fuel  change 
with  fuel  burnup  from  the  BOL  to  EOL.  The  changed  parameters 
govern  the  worth  of  control  rods  and  the  importance  of  any 
reactivity  insertion  or  removal.  The  transient  behaviour  and  the 
inherent  safety  features  of  the  reactor  are  also  dependent  on  the 
calculated  kinetic  parameters.  Hence,  the  reactor  will  have  differ¬ 
ent  behaviour  at  EOL  than  that  at  BOL.  Also,  the  neutron  flux 
throughout  the  reactor  increases  at  EOL  which  will  reduce  the  time 
of  isotope  production. 
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